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1 Introduction 

A tree is an acyclic structure where each node has zero to multiple descendent nodes and one parent node. A 
recombination tree is a special tree structure of which the size grows linearly. Therefore the investing decisions, 
if computed recursively, has time complexit}0 at most 0(n^), which is much more efficient than a general sim- 
ulation method which may cost exponential amount of time. For example, a recombination tree can help us to 
efficiently determine the price and the buy/sell timing for an American style option by comparing the deriva- 
tives value at each tree node with its children nodes. However, designing a recombination tree algorithm for 
modeling interest rates is far from trivial. Here are two examples: 

Hull and White {Si provided a heuristic two-stage method for constructing an interest rate tree based on the 
extended- Vasicek short rate model. In the first stage the algorithm builds the framework of the tree, and in 
the second stage the algorithm calibrates the tree to the current interest-rate term structure. The algorithm is 
designed for a short-rate model; hence the tree cannot be adjusted or updated according to the markets. Also, 
their method cannot deal with stochastic mean-reverting parameters, and there is no guarantee that the tree is a 
recombination tree especially when the volatility term of the short-rate process is a decreasing function. There- 
fore, Hull-White's algorithm is not a good candidate. 

Black, Derman, and Toy El (BDT hereafter) also provided a recombination tree algorithm for short rates. Their 
tree is constructed recursively and calibrated to zero-coupon bond volatilities and current interest rate structure. 
Though the BDT tree guarantees a recombination structure, the tree is not designed for a general Ornstein- 
Uhlenbeck process. Therefore, the BDT tree is not a good candidate either. 

In light of the fact that no existing tree algorithms can guarantee the recombination property for general 
Ornstein-Uhlenbeck processes, we propose a new recombination trinominal tree algorithm for our stochastic- 
splines model (this is formalized below). The idea is to modify a standard trinominal tree (Exhibit[T] [Top Left]) 
by adding extra branches at each node. First, we denote the tree structure in black color the center path. Then, 
given a node V directly above (below) the center path, define {V) the set of nodes containing V and all the 
nodes down to (up to) the center path. Then we modify the tree according to the following rules: (i) the center 
path remains unchanged; (ii) given a node V above the center path, we connect the node to all the descendant 
(children) nodes stemming from {V)\ (iii) given a node V below the center path, we connect the node to all the 
descendant nodes stemming from all the nodes from {V). The modified tree structure is shown in Exhibit[T][Top 
Right]. We will use the names, spanning nodes and spanning branches, to identify those nodes not on the center 
path, and branches not emanating from a center node. 



For the terminology of computational complexity please see Cormen et al, (2). 
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Exhibit 1: [Top-Left] Original Recombination Trinominal Tree; [Top-Right] Modified Tree; ; [Bottom] Simph- 
fied Tree 

The crucial key of modifying a standard trinominal tree is that we can further simplify and still keep the tree 
structure by adding sibling branches. Sibling branches are one-way streets through which we can only move 
up or down at a given time epoch, but not in both directions. A spanning node above the center path can reach 
the center path and all the nodes in between only by moving down through the sibling branches; a spanning 
node below the center path can reach the center path and all the nodes in between only by moving up through 
the sibling branches. As a result, by adding the sibling branches, each node can reach all but one descendent 
nodes via its sibling branches. So, in the simplified tree structure, each spanning node will have only one time 
transition descendant branch. The final tree structure is shown in Exhibit[T](Bottom). The algorithm is given 
below. The proof of the correctness of the algorithm for simulation is given in Section 3. 



2 Algorithm 

Now we give a full description of the algorithm. Let gjj denote the node on the center path at time tj, and let 
g{tj, j) denote the value of node gjj. Therefore, if g : R x Z i-> R is represented as a function, then it indicates 
the value of the node. Let gjj+t and g{tj, j + k) denote the A:-th node above center node gjj and the value of 
gjj+k respectively. Similarly, let gjj-k and g(tj, j - k) denote the fc-th node below center node gj j and the value 
of gjj-k respectively. Moreover, if we use capital letter G{tj, w), then it represents a random variable of the tree 
value at time tj. To shorthand the notation, the expectation value G{tj+i) conditional on the position of G{tj) is 
written as 




(1) 



Define the conditional expectation at node gjj to be 




(2) 
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Since the volatility term in stochastic- splines model is assumed to be a deterministic function, the conditional 
variance is the same for all nodes at a given time, i.e. at time tj, the conditional variance is 

Vj^Yai{G(tj^i)\G(tj)). (3) 

The idea of the recombination algorithm is to construct the center path first including the node values and branch 
probabilities, then determine the values of the spanning nodes, the probabilities on the spanning branches, and 
the probabilities on the sibling branches. The details are given in Algorithm[T] However, the algorithm shows 
that each tree is designed for simulating one coefficient process; if we have coefficient, we will need to build 
A'^ trees altogether if coefficient processes are correlated. After constructing the coefficient recombination "for- 
est", we can simulate the interest rate curve efficiently. 



3 Verification 

First we look at the^rif part of the algorithm and some notations. The first stage of the algorithm follows the 
standard Hull- White methodology and provides the backbone of the tree. Let gjj denote the node on the center 
path at time tj, and let g{tj,j) denote the value of node gjj. Therefore, if g is represented as a function, then 
it indicate the value of the node. Let gjj+k and g{tj,j + k) denote the k-th node above center node gjj and the 
value of gjj+k respectively. Similarly, let gjj-k and g{tj, j - k) denote the A:-th node below center node gjj and 
the value of gjj-k respectively. Moreover, if we use capital letter G{tj), then it represents a random variable of 
the tree value at time tj. To shorthand the notation, the expected value G(fj+i) conditional on the position of 
G{tj) is written as 



e[G(0+i)|G(0)1. (4) 



Now we move to the second part of the algorithm. The tree branches besides the central path are called spanning 
branches and spanning nodes. The second stage of the algorithm adopts the ideas of the law of total expecta- 
tions and the law of total variances to assign the values and probabihties of spanning nodes and branches. The 
procedure is done recursively. Therefore we just need to look at the cases when k - —\ and k - \. Given a node 
gj+ij+2 spanning from node gjj+\ at time tj, we denote the conditional expectation and conditional variance at 
node gj,j+\ to be Mjj+i and Vjj+i respectively. Denote the probability p to be the probability moving down 
from gjj+i to gjj and (1 - /?) to be the probability moving through the spanning branch from gjj+i to gj+\j+2- 

We can recall the law of total expectation which states 

E[X]^J]p(Y^y)E[X\Y^y]. (5) 

If we let 

X^E[G(tj^i)\G(tj)^gjj^i], (6) 
and Y denotes the random variable such that 

y _ f 0' if moving to the spanning brance 
[ 1, otherwise 

then 

Mjj.i =E[Gitj,,)\G(tj)^gjj,i] 

= pE[G{tj) I G(tj) = gjj] + (1 - p)E[G(tj^uj + 2) I G{tj) = gjj^i] 

= pMjj + (l-p)g(tj^uj + 2), (8) 
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Algorithm 1 Recombination Tree 



Require: MAXLEVEL (Tree Size) 



1: {Stage One - Center Path] 
2: for y = 1 to MAXLEVEL do 



3: 
4: 



Ax / 



round 



l M{l„j) \ 




<— /i X Axj 
<— (h + 1) X Ax^- . 
<— (h - 1) X Axj 

(M(tj,j)-g{tjj)y- M(iij)-g{tj,j) 

6 6VJ. 2V3V, 



iM(t,,j)-gitj.j}f 

iM{rj.j)-g(rj,j)f 
6V} 



M(tj,j)-g{tjJ) 
2V3V, 



7: end for 



8: {Stage Two - Spanning Branches} 



for i = 1 to MAXLEVEL do 

for j = 1 to 2 X J - 1 do 

{Move to the j-th vertex below the center path} 

Find X e [M,j_i - ^/V'i, Mij^i + ^fVi] where 



X - Mij Vj + {Mij - Mij^i)^ -{x- Mij.r)2 



13: g{ti+uj-2)<^x 

15: p{gi,j-UgM,j-2) <^l- '^^^ 

16: {Move to the j-th vertex above the center path} 

17: Find X 6 [M,j+i - yfVi, Mij+\ + y[Vi\ where 

X - Mjj+i Vi-{x-Mi,j+if 

X - Mij ~ Vj + (Mij - Mij^yf - (x - M,- ,+i)2 



18: g(f,+i,; + 2) <- X 

19: P(gi,j+X,gi.j) ^ 

20: ffej+l>§f+lj+2) <- 1 - — 

21: end for 

22: end for 
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which shows 

P §{tj,uj + 2)-Mjj ■ 

The task of deriving the relationship between g{j + 1,7 + 2) and p from the law of total variance is more 
complicated. We will show the the result first, then break into each part. Recall the law of total variance which 
states 

Var(X) = E [Var(X | Y)] + Var(E[X | Y]) . (10) 

Similarly we let 

X^E[G(tj^i)\G(tj)], (11) 
and Y denotes the random variable such that 



0, if moving to spanning brance 

1 , otherwise 



(12) 



If the following statement is true: 

y^r(E[G(0+i) I G{tj) = gjj^i]) = Efy^r(EfG(0+i) I G{tj) = gjj^i] | y)] 

+ V^r[E[E[G(tj^i) I G(tj) = gjj^i] | y]) 

= pVj + p{Mjj - Mjj^i)- + (1 - p)ig(j, ;■+!)- Mjj^i)-, (13) 

then we have 

P Vj + {Mj,j-Mj,j,,)^-(g(jJ+l)-Mjj,,y~- 
Examining the first term, pVj, on the right-hand-side of equation ( fTSl ). 

V^r (e [Gitj^O I G(tj) = gjj^i] I y = O) = y^r {g(j + 1 , ; + 2)) = (15) 
since there is only one choice moving from gjj+\ to g j+\j+2- On the other hand, 

V^r (e [G(0+i) I G{tj) = gjj^i] \Y^l)^Vj, (16) 
and we know this value recursively. So 

E[y^r(E[G(0+i) I G{tj) = gjj^i] I y)] = (1 - p)V^r{E[G(tj^i) | G(tj) = gjj^i] | 7 = o) 

+ pV^r{E[G(tj^i) I G(0) - I y = l) 

= /7y^r(E[G(0+i) I G(0) = gjj+i] I i- = l) 
= pV^i- (17) 

Next, the second and third terms. Since 

e[e[e[G(0+i) I G(tj) = G{tj,j + 1)] I y]] = e[G(0+i) I G(tj) = G(tj,j + 1)] = Mjj^u (18) 

we have 

y^r(E[E[G(0+i) I G(tj) = gjj^i] I y]) 

= (e[e [G(0.i) I G(tj) = I y - l] - Mjj,,f 

+ (l-p) (e [e [G(0.i) I G(tj) = I F = O] - Mj,j,,f 

= (e [G(ry+i) I G(tj) = ^,;,] - Mjj^if + (l-p) [gU + 1, + 2) - M,;,.+i)' 

- p (M;.,. - M;.,.+i)' + (l-p) {g{j + 1, + 2) - My,,.+i)' . (19) 



5 



Now we have two equations and two unknown p and g(tj+i,j + 2) in the following 

C _ g(tj+i,j+2)-Mjj^i 

J ^ g(tj^,.j+2)-Mjj 

1 ^ V/-(gO'+lj'+2)-A^;.,>lf ■ 

However, p must be a number between and 1 . And we now show that the equations indeed yield a solution 
such that p e [0, 1]. First, the case where Mjj+i > Mjj and write 

/iW = TJ^ (21) 



and 



Vj + (Mjj-Mj,j,0^-{x-Mjj,ir 



^(^)- .^ .^.^ ^2_.._......2 - (22) 



Since My^+i > M^^-, for any x > Mj^ 



X - Mj ;+i 

< — < 1 (23) 



and /i(x) is continuous and monotonically increasing. On the other hand, for any x e [Mjj+i, Mjj+i + ^fVj], 

< ^— < 1 (24) 

Vj + (Mjj-Mjj,,)^-{x-Mjj,i)^ 

and f2{x) is a continuous and monotonically decreasing function. Since 

MMjj^i) = 0, MMjj^i + ^/v'j)>Q (25) 

and 

fiiMjj^i) > 0, MMjj^i + V^;) = 0, (26) 

we know that there must exist a unique x € [My.j+i, Mjj+i + ^JVj] such that 

/i(^) = /.(^) = pe[0,l]. (27) 

Alternatively, the proof is similar for the case when Mjj+i < Mjj except the solution exists in [Mjj+i - 
yfVj, Mjj+i]. The uniqueness and existence of the solution p and g{tj+\,j + 2) help us solve the equations 
fast. 

4 Application: AA Rated Callable Corporate Bonds 

A callable bond is a bullet bond with an embedded American-style issuer call option. The option gives the 
issuer the right to refinance (when interest rates are low, for example), therefore the price of a callable bond 
is lower than the price of a bullet bond with the same coupon rate and maturity. To price the option, we can 
assume the underlying bullet bond price process follows a log-normal process. However, the method does not 
work if we simply want to price the AA corporate callable bond with its option feature. A better solution is to 
model the bond yield term structure instead of the individual bond price. Further discussions can be found in ||4| . 

The purpose of this application is to demonstrate how to simulate AA bond yield curve by the recombination 
trees constructed according to the algorithm described in previous sections. Particularly we consider only AA 
corporate bond universe, and model the curve under the stochastic-splines-model framework in Q. The model 



6 



is constructed by assuming that the coefficient processes follow the exponential Ornstein-Uhlenbeck processes 
under the real-world measure, which later can be transformed to a risk-adjusted measure by imposing proper 
no-arbitrage conditions. Once the trees are built under both measures, we can simulate the bond yield curve 
under the real-world measure to price bullet corporate bonds, and simulate the bond yield curve under a risk- 
adjusted measure to price bond embedded options. Then callable bond price will be the bullet bond price minus 
the option value. 

In Exhibit in we demonstrate those recombination tree results under the real-world measure. In Exhibit [3] 
we demonstrate those recombination tree results under a risk-adjusted measure. The black branches represent 
the center paths; the blue branches represent the spanning branches; the green branches represent the sibling 
branches that can only moves down; the red branches represent the sibhng branches that can only moves up. 
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Exhibit 2: Real-World Recombination Tree Results based on 08/31/2010 AA Corporate Bond Yield Curve 
Dynamics 
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Exhibit 3: Risk- Adjusted Recombination Tree Results based on 08/31/2010 AA Corporate Bond Yield Curve 
Dynamics 
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